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This work analyses the temporal coherence and correlation of counterpropagating twin photons 
generated in a quasi-phase matched nonlinear cristal by spontaneous parametric-down conversion. 

We find out different pictures depending on the pump pulse duration relative to two character¬ 
istic temporal scales, determined respectively by the temporal separation between the counter¬ 
propagating and the co-propagating wavepackets. When the pump duration is intermediate between 
the two scales, we show a transition from a highly entangled state to an almost separable state, with 
strongly asymmetric spectral properties of the photons. 
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Introduction 


Spontaneous parametric down-conversion (SPDC) occurring in x ^ media is one of the most accessible sources 
both of entangled photon pairs and of single photons, heralded by detection of the partner. The microscopic process, 
where a high energy photon of the pump laser splits into two lower energy photons, is ruled by conservation laws 
(energy, momentum, angular momentum, polarization), which are at the origin of a wide range of quantum correlations 
between the members of the pair. 

In standard co-propagating configurations, the two-photon state is characterized by a high dimensional entangle¬ 
ment because a quantum correlation is present over huge temporal and angular bandwidths. The temporal correlation 
was historically the first one to be studied pQ: a down-conversion event can take place anywhere along the crystal, 
so that the arrival time of the twins is not known. However, the members of a pair, generated at the same point, 
propagate nearly in the same direction, and exit the crystal almost simoultaneously A small uncertainty in their 
temporal separation is present because of their different group velocities (type II) or because of the group velocity 
dipersion (typel), and can be reduced to its smallest limit (the optical cycle) when the spatial degrees of freedom 
are properly controlled 0 , 0 - Such a short correlation time results in a high-dimensional temporal entanglement 
[5]. Its spectral counterpart is the huge spectral bandwidth of SPCD emission, and the high dimensional spectral 
entanglement of SPDC photons[fL High-dimensional entanglement offers relevant opportunities in view of broadband 
quantum communication schemes, but can also be regarded as a negative feature, because it affects the purity of 
heralded single photons. 

This work considers a non-conventional configuration, where one of the down-converted photons is generated in the 
backward direction with respect to the pump laser, in a periodically poled crystal (Fig0- Although predicted almost 
fifty years ago [7 , counter-propagating down-conversion has been only recently realized 0i, thanks to technical 
advancements in achieving the sub-micrometer poling periods necessary to phase-match the interaction m ■ 
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FIG. 1: (Color online) Geometry of the counterpropagating down-conversion (see text). 


Counterpropagating PDC presents unique features, as the presence of a threshold for the pump intensity beyond 
which coherent parametric oscillations take place m, thereby the name Mirrorless Optical parametric Oscillator 
(M0P0)[8]. In a related work [12], we study the quantum correlation of counter-propagating twin beams close to 
the threshold; here, instead, we focus on the regime of spontaneous photon pairs production, well below threshold, 
and analyse the temporal quantum properties of counterpropagating twin photons generated in a purely collinear 
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configuration. 

A second peculiar feature of the MOPO is the narrow spectral bandwidth of emission (the backward propagating 
wave can be even more monochromatic than the pump laser [8]). In the quantum domain, as pointed out in m, 
counterpropagating SPDC can generate highly monochromatic photon pairs in an almost separable state, which makes 
it a promising source of high-purity heralded single photons. 

In this work we provide a detailed theoretical analysis of the effects of the spectral properties of the pump laser on 
the degree of entanglement of the state, identify the physical conditions under which the state may become separable, 
and provide a consistent interpretation of the transition from an entangled to a separable state. 

In particular we will show that the system dynamics is governed by two well separated time scales: a long one , 
related to the temporal separation of counterpropagating waves, which is on the order of the transit time of light 
along the crystal (tens of picoseconds), and a short one related to the temporal separation of co-propagating waves 
due to their different group velocities (order lps or smaller). When the pump pulse duration is intermediate between 
the two scales, we will show that the two-photon state becomes separable, and remains separable for a wide range 
of pump durations, whilst it has a high degree of entanglement in the two opposite limits. Notice that such a 
difference of time scales occurs naturally in the counterpropagating configuration, for basically any kind of material 
and tuning condition. This is quite different from the co-propagating case where separability of the state requires 
special operational points [26]. 

In addition, we shall investigate the coherence properties of the SPDC photons taken individually, showing a 
transition from a symmetric state, for a long pump pulse, to a highly asymmetric state for a short pump pulse. In 
particular, in the regime where the state is separable, the spectrum of the signal turns out to reproduce the spectrum 
of the co-propagating pump laser, while that of the backward propagating idler is entirely determined by the crystal 
properties. 

The paper is organized as follows: Section [I] introduces a quantum model for counterpropagating PDC; Sec. [I I 
characterizes the quantum correlation of twin photons in the spectral domain, while seeing is devoted to their 
coherence properties. SecJTV| quantifies the degree of entanglement of the state via the Schmidt number. Finally 
SecfV| gives an interpretation of the transition from entanglement to separability by analysing the correlation in the 
temporal domain. 


I. THE MODEL 

The starting point of our analysis are the equations describing the propagation of three interacting pump, signal 
and idler waves along a slab of periodically poled second order nonlinear crystal (Fig.l). We shall consider here only 
collinear propagation, either assuming that light is collected only at small propagation angles with respect to the 
pump, or because of a waveguiding configuration. In a crystal with periodic inversion of the nonlinear susceptibility, 
the momentum conservation in the three-wave interaction is replaced by a less restrictive conservation law (quasi¬ 
phase matching) [I?! , which includes also the momenta of the nonlinear lattice, where A is the poling period (for 
example m = ±1, ±3... for a simple poling). Since the effective nonlinearity is higher for lower orders, one usually tries 
to phase-match the first order m = ±1 interaction. The counter-propagating configuration, in which one wave (say 
the idler) is generated in the backward direction with respect to the pump laser (Fig.l), thus requires a poling period 
on the order of the pump wavelength, because the pump momentum needs to be almost entirely compensated by the 
grating momentum ko = In these conditions, the central frequencies cj s , uji = uj p — uj s of the down-converted 
wavepackets are determined by quasi-phase matching at the central pump frequency uj p 

27r 

koi = kQp ^ (1) 

where k G j = /c, j = are the wave numbers at the three central frequencies. We shall mostly focus on 

the commonly realized type O interaction [8], where the three waves have the same polarization, but we leave the 
formalism quite general. Hence the subscript j in the wave number may refer to dispersion relations for either the 
ordinary or extraordinary wave, including thus type II or I PDC. 

We than introduce the positive frequency part of field operators (with dimension of photon destruction operators) 
for the three wavepackets as 


A s (Cl, z) 

= e +ik A Q >a s (Ct,z), 

(2a) 

Ai(Q,z) 

= e- ik ^ n)z di(n,z), 

(2b) 

Ap(yi, z) 

= e +ik A^ z d p {n,z), 

(2c) 
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where capital Q is a frequency offset from the carrier frequencies, and kj(Q) are the wave numbers at frequency Uj + Q. 
In this definition, the factors e ±lkj ^ z account for all the effects of the linear propagation along the medium. Hence 
the operators dj have a slow variation along the crystal because they evolve only under the effects of the nonlinear 
interaction. Their coupled equations of propagation can thus be written as ( see [12] for a more detailed analysis) 


-^-a s (n,z) = x j dfl'a p (fl + z)a\(Q !, z)e lV (^ l ^) z ^ 

(3a) 

z) = —x J dn'a p {n + n / ,z)al{Cl',z)e- w(Q '’ Q)z , 

(3b) 

Ta p ( 0 , z) = -xj z)di{n - n', Z ) e «>(n,n-n')z ; 

(3c) 

where x is proportional to the effective second order susceptibility of the crystal, and only the first order terms ±1 in 
the Fourier espansion of the periodic nonlinear suceptibility have been retained (namely order -1 for signal and idler, 
order +1 for the pump) . In these equations 

O 7 r 

v(n, n') = k s (n) - ki(tf) - k p (n + n')' + — 

(4) 


is the effective phase mismatch that rules the efficiency of each elementary down-conversion process, where a signal 
and an idler photon of frequencies uj s + £ 2 , uji + Q' are generated out of a pump photon of frequency uj p + Q + Q f . 
Notice that, a part from the different form of the phase matching Q , the only formal difference with the usual 
co-propagating case (see e.g. m) is the minus sign appearing at r.h.s of (|3|b for the counterpropagating idler. As 
we shall see in the following, however, this minus sign leads to very relevant physical differences. 


A. Low-gain limit 

In a parent work m we analyse these equations in a generic gain regime, including also the region close to the 
MOPO threshold. In this work we instead focus on the ultra-low gain regime, much below the MOPO threshold, 
where photons pairs are generated by purely spontaneous down-conversion. In this regime, the depletion of the pump 
beam can be neglected and the pump approximated by a constant c-number field, corresponding to the pump pulse 
at the crystal input face 


d p (Q, z) —>• a p (fi, z) ~ cq,(Q, z = 0) (5) 

The strength of the parametric coupling is then described by the dimensionless gain parameter 

g = V2nxa!p(t = 0)l c , (6) 

where a p (t = 0) is the peak value of the pump temporal profile. Notice that in the limit of a monochromatic pump 
mm g ss 7r/2 represents the threshold for the MOPO. Conversely, in the limit g <C 1 Eqs. ©> can be solved 
perturbatively. Namely, we write the formal solution of (|3|, starting from the boundary conditions: 

d s (Q,z = 0) = (7) 

d i(i},z = l c ) = a- n (Q) (8) 

determined by the input signal and idler fields, entering the crystal from the left face of at z = 0 and from the right 
face at z = Z c , respectively (Fig.l) . By solving iteratively, a perturbative seriers of powers of g is obtained. By 
keeping only the first order terms in g 1 , one obtains a Boguliobov linear trasformation that links the output to 
the input operators: 


d,s (^s 5 % - 

lc) 

(9a) 

Ws) + 

J dfliip(£l s , f2j)a‘ nt (Oj) , 

(9b) 

d{ (fos •> % = 

0) 

(9c) 

d- n (Oi) + 


(9d) 
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with the biphoton amplitude given by 


5 ^z) — 


X 


nr —cq,(u s H~ £4*) 

V Z7T 

V(Q s ,Qi)l c 


. X>( 12 s Mi)lc 


where 


( 10 ) 


dp(U) 



dt c iQt <*(t) 
\[2iT a p {t = 0) 


( 11 ) 


is the Fourier profile of the pump pulse at the crystal input face, normalized to its temporal peak value. Notice that 
Eqs.(|9| define a unitary tranformation only up to first order in g. In the following, the input signal and idler field at 
the left and right end faces of the crystal will be taken in the vacuum state. 

It is worth remarking that the quantum field formalism here employed can be replaced by an equivalent state 
formalism (see e.g [4] for more details) where the state evolves instead of the field operators. By applying the 
trasformation 0 to the input vacuum state, one obtains at the output the usual state 

10 )out = |0) + I J dQ. s dn i ^(Q. s ,Q, i )al(Q. s )al(n i ) |0) , (12) 


describing the superposition of the vacuum state |0) and of a two-photon state, where the photon pair can be generated 
in any of the Fourier modes with probability amplitude In this respect, the formalism used here is 

equivalent to the one employed in mm- 


II. SPECTRAL BIPHOTON CORRELATION 


This section is devoted to the analysis of the biphotonic correlation in the spectral domain. Precisely, we focus on 
the probability amplitude (A° ut (D s )A° ut (f^)) of finding a pair of photons at frequencies at the crystal output 

faces. Using the input-output relations ([9| and the definitions ©: 

(A° ut (^)^> ut (^)> = e ifc «( n *)‘^(n„ no. (13) 


with ^ given by Eq.(10). As usual, the biphoton correlation is the product of two terms: i) the pump spectral 
amplitude a p (fl s + U7)7 reflecting the energy conservation in the microscopic process, and ii) the phase matching 
function smc(Vl c /2)e ~ lT>lc ' 2 , reflecting the generalized momentum conservation. Concerning the latter, we can expand 
V(Q s ,Qi) in Eq.Q in power series of the frequency shifts from the carriers. Down-conversion spectra are typically 
narrow M, as will become also clear in the following, so that one is allowed to retain only terms up to first order 


v(n s ,Qi)j 


w- 


tos 


Q 


gvm 


k' p )n s -(k' + k' p ) fq 


Hi 

Hi 


gvs 


(14) 

(15) 


where the zero order term vanishes because of Eq. 0 - and k[ = dkj / d(jj\^ =u , j = s,i,p. We thus see the appearence 
of the two characteristic temporal scales: 


Tgvm 


:=a 


gvm 


2 Hgfp Vg S 


•= Q' 

^ L gvs 


-1 



(16) 

(17) 


where v 9 i = 1 /k[ are the group velocities of the three wavepackets at the central frequencies. The first scale 
[Eq. (pU|) ] describes the ” small” temporal separation between the co-propagating waves due to their group velocity 
mismatch (GVM) . The second one [Eq. <EZ1>] accounts for the ” large” temporal separation of the counter-propagating 
pump and idler waves, which is ruled by the time needed by the pulse centers to cross the crystal. Closely related, 


T gvs 


= a 


-i 


gvs 


1 

2 


'lc_ 

, v gs 


lc_~ 

v gi . 


(18) 
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describes the characteristic temporal separation between the arrival times of an idler and a signal photon at their exit 
faces. Clearly, since group velocities are close, r gvs ~ Tg VS , while r gvm <C Tg VS ,r gvs , and 


T gV m 

V = -r- 

r gvs 


Q' 

^ ^gvs 

f7gvm 


< 1. 


(19) 


Therefore, the phase matching has two well separated scales of variation: as a function of the signal frequency it decays 
on the broad bandwidth f2 gvm , while as a function of the idler frequency it decays on the narrow bandwidth fl gvs . 
Plots of the parameter rj , for periodically poled KTP (potassium titanyl phosphate) and LiNbOs (lithium niobate), 
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FIG. 2: (Color online) Ratio r] — r gvm /rg VS for periodically poled KTP and LiNbOs, pumped in the infrared or visible, for 
type 0 e —y ee down-conversion. Point A is KTP pumped at X p — 821nm, with A po i = 800nm, A s = 1141nm, A i = 2932nm, 
corresponding to the experiment in [8]; B is the zero GVM point for the KTP at X p = 821nm, corresponding to A po i = 290nm, 
A s — 3523nm, A* = 1071nm. C is a LiNbOs slab pumped at X p = 527.5nm, for degenerate PDC at A s = A i = 1055nm, with 
A P oi = 236nm. 


are shown in Figj2j where A,B,C are the points that will be used as examples in the following. 

Finally, a third relevant scale is the pump spectral bandwidth. For a coherent Gaussian pump a p (t) = 
a p (0) exp —^ 2 , the pulse duration r p is the inverse of the bandwidth 


t d = 


1 


( 20 ) 


Depending on the pump bandwidth relative to the spectral scales of phase matching, different physical situations 
arise. The three relevant possibilities, depicted in Figj3j will be studied separately in the following. 


i) Limit of a CW pump: 

We assume a narrowband pump pulse, such that 

Tp Tg Vg T gvm , Or AOp ^ gV s ^ ^gvm • (21) 

This limit corresponds to a pump pulse that in the £ direction is much longer than the crystal slab, and for a crystal of 
some mm length requires a pulse duration of hundreds of picoseconds or longer. In this limit the pump spectral profile 
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FIG. 3: (Color online) Biphoton correlation \'ip\ [Eg. ( |10| )] in the plane (TEfU) , in various pumping regimes. Example of a 
4mm PPKTP, pumped at 821.4nm, corresponding to the point A in Eigf2j with Tg VS = 25.2 ps, r gvm = 0.27 ps. (a) Quasi CW 
pump pulse r p = 253ps. (b) Intermediate pump pulse r p = 1.1 ps. (c) Ultrashort pump r p = 0.03 ps . Note the different scales 
of the plots 10 11 —> 10 13 Hz. 


a p (fl s + Qi) is much narrower than the phase matching bandwidths, and the geometry of the correlation is dominated 
by energy conservation, which requires that the twins are generated at symmetric frequencies Q s + = Q p « 0. As a 

consequence, the biphoton correlati on (|I()| ) has a sharp maximum along the diagonal Q s = —Q*, as shown by Fig(3^i. 
Indeed, as derived in the Appendix |A1[ in this limit the correlation is well approximated by : 

ip(Q. s , Q*) ~ —|=d p (fi s + f2j)sinc ( 4F- ) 

V 27r \' ^gvs J 

- -jL=a p (tt s + Oi)sinc ( ) e ^ 

v 27 r \'^gvs / 

ii) Limit of an ultrashort pump pulse: 

We consider here the limit: 

Tp <C Tg vm , T gvs , or AQp Qgvm? ^gvs ? (24) 

where the pump pulse is not only shorter than the crystal length, but also shorter than the average separation between 
the pump and signal wavepackets due their GVM. In our examples this corresponds to duration shorter than 100 fs. 
In these conditions, the pump spectral profile a p (fl s + f^) decays slowly with respect to sincD(Q s , f^)Z c /2, so that 
the geometry of the biphoton correlation is dominated by the phase matching, i.e. by the momentum conservation. 
As a result (Fig)^) the biphoton correlation takes the approximated form 


( 22 ) 

(23) 


— r -[^s(l ^ 7 )] sine 

V 2tt 


to* 


+ 


Q ? 


O O' 

i igvm * ^gvs 


fi gvm ^gvs 


5 


V2tt 


-a 


?7 


sine 


ft* 


+ 
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When plotted in the plane (ft, . Q s ). the function shows a sharp maximum along the line 

n a 






L _ gvm 

n' 


(25) 

(26) 


(27) 


where phase matching occurs (see Eq. (p~5])) , 


and very asymmetric spectral properties of the signal -idler photons. 


iii) Intermediate pump pulse: 

The intermediate case, where 


'Tgvs ^ Fp ^ T”gvrm Qg VS AQ p f^gvm 5 


( 28 ) 
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is the most peculiar one, because the biphoton correlation may approach a separable function of (FigJS^). 

First of all, we remark that the limit ( [28] ) is strictly realized only for r\ = T gvm /Tg VS 0, i.e for a vanishing group 
velocity mismatch between the pump and the signal. This condition is favourable to separabilty, because as 77 0 

the phase matching function tends to become a stripe parallel to the axis (see Eq. (27) ), but it is not a sufficient 
one, because of the role of the pump profile in Eq.([l0|). H owev er, provided that the the pump spectrum satisfies the 
intermediate limit (28), it can be demonstrated (Appendix |A 1 | ) that the the biphoton amplitude ( 10 ) approaches the 
factorized form: 


xjj(Q s ,Qi) —» _ a p (fl s )e ^s vm x sine 

v 27T 



Pi 

1 O 7 

e s vs , 


(29) 


i.e. it becomes the product of a function of reproducing the pump profile, and a function of f^, corresponding to 
the phase matching profile. This describes a nonentangled biphoton state, with the signal photon generated in the 
same spectro- temporal mode as the pump, while the spectral mode of the idler is dictated by the phase matching 
”sinc” function of width f 2 gvs . 

This qualitative picture will be confirmed by the evaluation of the Schmidt number in Sec. [TV| a nd will be further 
interpreted and discussed in the light of the temporal correlation of biphotons described in Sec|V| 


III. SPECTRAL COHERENCE OF COUNTERPROPAGATING PHOTONS 


This section is devoted to the marginal statistics of individual signal and idler photons. The focus is on their 
spectral coherence properties, studied by means of the first order coherence functions 


GW(f2,fy) = 

Gj 1 } (n,n') = {Al out (n)A° ut (n')) 


(30) 


(where a propagation phase factor is present in the first definition just for convenience of notation). From the 
input-output relations © one has: 


G«(ft s ,ft' s ) = J d^v^WOW), 


(31) 

(32) 


i.e. the coherence functions are convolution integrals over the biphoton amplitude given by Eq(fIo|) . The knowledge 
of the is sufficient to determine all the statistical properties of the marginal distributions. For example, the 
autocorrelation of the light intensities Ij = A-A, is given by 


(/,(«)/,(«')) = s(q - + (ijimW)) + A(^) 


(33) 


where = G^(Q, H). This relation, which is a consequence of the factorization theorem of Gaussian moments, 

is typical of thermal-like statistics. As a matter of fact, the marginal distributions of the output signal-idler light 
are thermal-like Gaussian, when there is vacuum at the input. In the low-gain regime considered here, the dominant 
term is the first one, i.e. the ” shot-noise” term ^-correlated in frequencies. Therefore, as well known in this regime 
the statistics of photon counts in each arm is Poissonian. 


On the other side, the convolution integrals in Eqs. (31), (32) indicate that an autocorrelation of spectral fluctuations 
inside each individual signal or idler wave exists because of second order processes, that involve the probability 
amplitudes of generating at two pairs of photons. 

In the following we shall illustrate the three relevant cases. The coherence functions will be evaluated both nu¬ 
merically (F igj4| and analytically. In the first case, the complete Sellmeier relations [18] will be used to compute the 
integrals in ( |31j ), (32), while the linear approximation for phase matching will be exploited to derive approximated 
analytical formulas. 

i) Limit of a CW pump: 

Column (a) of Fig|4] shows an example of the signal and idler coherence functions in the plane (fi, 12'), numerically 
computed in the case of a long pump pulse r v ~ 14rg VS . 
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FIG. 4: (Color online) The coherence functions |G(Q,fT)| of the forward signal and backward propagating idler are plotted 
in the upper and lower row, respectively, for different pumping regime . Column (a) Quasi CW pump pulse r p = 353ps. (b) 
Intermediate pump pulse r p = 1.1 ps. (c) Ultrashort pump r p = 0.03 ps. Same KTP crystal slab as in Fig|3] (point A in Fig. 
[5J, with Tg VS = 25.2 ps, Tgvm = 0.27 ps. Note the different scales in the panels 


In the limit r p Tg VS , approximated expressions for the coherence functions canj3e calculated by inserting the 


formulas for the biphoton correlation (22) and (23), valid in this limit, into Eqs. (31) and (32), respectively, and 
performing the simple integrals. After some passages: 


gF(C,0') (34) 

-Tp»^v S 2^, _ ^ 5 2 sinc 2 (35) 

where T p (£t) = f ^e lQt \a p (t) / a p (0)\ 2 is the Fourier transform of the pump intensity profile. These approximated 
formulas have been checked with the numerical results and show an excellent match. These results may be considered 
the more refined version of the much simpler CW model analysed in mm, with the narrow peak i p (Q' — Q) being 
the finite counterpart of the singular Dirac S appearing in the strictly CW pump model [12] . 

For a quasi-CW pump the counter-propagating signal and idler photons are predicted to have identical spectral 
coherence properties. In particular, by looking at the functions along the diagonal Q' = Q we see that their 
spectra (Jj(Q)) = (A-(Q)Aj(Q)) 


(i s (^)) = </i(n)) y ^ sinc2 (fT - ) ’ ^ 

are identical and entirely determined by the narrow bandwidth of phase matching Q gvs . This bandwidth is in turn 
the inverse of the characteristic separation r gvs between the arrival time of an idler and a signal photon at their 
crystal end faces, which rouhgly corresponds to the long transit time of light along the crystal slab, because they 
propagate in opposite directions. As already noticed in m, and as will be further discussed in Secjvjthis is clearly 
a big difference with the copropagating case. There, the temporal uncertainty between the arrival times of the idler 
and signal photon is short, because determined at most by the group velocity dispersion or mismatch, which results 
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in the huge down-conversion bandwidths that characterize the standard co-propagating configuration. 

On the other side, when studied as a function of Q f — Q the G T) gives the characteristic size of spectral fluctuations, 
i.e. the spectral coherence length. This is determined by the pump bandwidth, more precisely by the width \/2Afl p 
of X p (fy — fi), which is much narrower than the spectral bandwidths f2 gvs . We can heuristically estimate the number 
of modes by counting the number of coherence length contained in the spectrum: therefore, for such a long pulse we 
expect each signal and idler photon to be generated in a highly incoherent and multimode state, with the number of 
modes oc = -^ E ~. 

ZAiir) Tcr v s 


ii) Ultrashort pump pulse: 

When the pump pulse shorten below the transit time r^ vs along the crystal slab, the spectral properties of the 
counterpropagating idler and signal change drastically, becoming strongly asymmetric. First we consider the case of 
an ultrashort pulse, r p T gvm (i.e. such that pump and the signal tend to split apart during propagation). The 
asymmetry between the forward and backward propagating photons can be clearly appreciated in the third column 
of Fig. [4j which plots their coherence functions for r p ~ 0.1r gvm . 

Approximated expressions for the coherence functions are derived also in this case, by using the limit behaviour of 
the biphoton correlation described by Eqs. (25) and ( [26] ) . With some calculations: 


GW(9,9') Tp ^ 


n 2n/ 

T’p'CT’gvm O “ "gVS 


<5 p [f2(l — r/)]\ 2 sine 


Cl' -Q\ 


Q„ 


• / n' -n \ 

\ ^gvm J 


(37) 


This formula predicts that the spectrum of the forward propagating signal 

9 2 K 


{i*m = 


i« P [^(i-»?)]r 


(38) 


is a replica of the pump spectrum with a scale factor on the order unity. The coherence length of 

the signal (the characteristic size of spectral fluctuations) is instead determined by the width of the narrower sine 
function, l co h,s ~ ^ gV m- From this picture we thus expect that the signal photon, when detected independently from 
its twin, is in a incoherent multimode state, with the number of modes oc — • 

In a similar way, for the idler photon we get: 




Tp<T g , 


g 2 n 


gvm 


Oi p [—12 


1 — rj. 


x sme 


"Tf 


gvs 


(39) 


This formula predicts an idler bandwidth much narrower than the pump, precisely it predicts that the idler spectrum 

k' —k' 

follows the pump spectrum with a scale factor = k'+k' ^ 1 • The coherence length of the idler is l co h,i ~ ^ gvs > 
so that the number of temporal modes is predicted to scale as (1 — , which is the same number as 

V 1 gvs V 1 V)^ £ gvm 

for the signal (as it must be because the signal and idler are the two members of the same entangled state, and their 
reduced states must exhibit the same Schmidt dimensionality, see next section) . 

Notice that this particular scaling of the bandwidths of the forward and backward propagating waves with the pump 

/c ^ j ( 

bandwidth is well known in the literature concerning the MOPO. There, the same scaling factors, k T +k ? for the 

k' —k' 

forward-propagating signal and fe ? +fe , a for the backward propagating idler, are predicted to occurr B®, by using 


arguments based on the phase-matching characteristic of the process. Here, however, the analysis concerns the 
quantum properties of the single photons generated well below the MOPO threshold. Moreover, at difference with 
the classical analysis in [8], such a scaling with the pump spectrum is predicted only in rather extreme conditions, 
corresponding to an ultrashort pump pulse r p <C T gvm . Notice that this limit imposes a precise and not trivial 
constraint on the minimum observable bandwidth of the idler photon: the behaviour described by Eq.(39) is indeed 
realized only for r p <C r gvm , or for >> ^ gV m, so that it requires that the idler bandwidth 




Tj Tj 

- AO p — 9g V m ftgvs 
1 — rj 1 — 7] 


(40) 


iii) Intermediate pump pulse: 

When r gvm « r p « Tg VS , the properties of the twin photons are actually intermediate between the two former cases, 
with the forward propagating signal photon replicating the pump spectrum, while the coherence properties of the 
backward propagating idler are determined by phase matching. These features are clearly exhibited by the central 
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column (b) of Figlil which plots a numerically computed example of the coherence functions for r p = 0.04r' « 4 77 


1 gvm) 


« - - t' gVD 

o the transit time along the slab, but long enough that GVM does not play a relevant role. 

The observed features are a straightforward consequence of the separable form (29) of the biphoton amplitude which 
holds in this limit. Indeed, by using Eq.(29), in the limit r p /r^ vs —> 0 ,T gvm /T p 0 we obtain: 


GW(n,n') 

G?\n,n') 


o 2 Q 

y * L gvs 


a* [ft (1 - T})] a p [O' (1 - T})] 


> 9 2 t p • 

-» ^-L sinc 


V2n 


ft 


a 


gvs 


sine 


ft' 
ft~ 


gvs 


(41) 

(42) 


Thus in this case the signal spectrum is a replica of the broad pump spectrum I s (ft) oc \a p [12 (1 — 77 )] | 2 , while the idler 
spectrum is determined by the much narrower phase-matching function Ii(ft) oc sine 2 (fp— Precisely, the signal 

spectrum is described by the same formula (38) as in the ultrashort pump case, while the idler spectral properties 
are described by the same formula (36) that holds in the CW pump limit. However, notice that in the present case 
the coherence properties are remarkably different, as the two coherence functions are perfectly symmetrical along the 
two diagonals ft ± ft': as can be easily inferred from Eqs. (41) and (42) the two coherence lengths are l co h,s ~ A12 p 
and l co h,i ~ ^gvs 5 he. they are equal to the respective spectral widths. This is in accordance with the separability of 
the biphoton state, which corresponds to single-mode, almost coherent reduced states for each of the two twin photon 
taken separately. 

We conclude this section observing that the results (36), (40) and (42) implies that in any pumping regime the idler 
bandwidth cannot be narrower than the phase matching bandwidth U gvs , a limitation that arises from the imperfect 
momentum conservation due to the finite length of the crystal slab. 


IV. SCHMIDT NUMBER OF ENTANGLEMENT 

So far our considerations about the number of modes and the degree of entanglement of the system have been 
qualitative. A quantitative measure of the entanglement is offered by the so-called Schmidt number mm, which 
is recognized to give an estimate of the number of Schmidt modes partecipating in the entangled state, i.e. of the 
effective dimensionality of the entanglement ED- First of all, as usual, we consider the state conditioned to a photon 
count 


\4> c ) = ■'{ dn s dniii>(n s ,fy)at(n s )al(fy)\o) 


(43) 


where with respect to the true output state ( 12 ), the vacuum term has been dropped. Then, we introduce the Schmidt 
number, as the inverse of the purity of the state of each separate subsystem 

1 1 


JC = 


T r{pl) Tr{pf} 


(44) 


where p s , pi are the reduced density matrix of the signal and idler , e.g. p s = Tri{|0c) {(j>c |}- For a two-particle state 
of the form (43), the Schmidt number can be calculated via an integral formula, as e.g derived in [4] (see also 0 ) , 


K = 


B 


where 


V = y dflG^\n,fl) = J dQGpiQ.Q) 

j dn J dn' G^in ,^) 2 = J dn J dn' G) 1 } (fi,fi') 


b = 


(45) 

(46) 

(47) 


As can be easily checked, J\f is the espectation value (first order moment) of the photon number operator Nj = 
f dftlj(ft) in either the signal or idler arm 


M = (N a ) = (Ni) 


(48) 
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The quantity at denominator is instead linked to the second order moment of the photon number. By performing the 
integral of Eq.(33) over the two spectral arguments, one gets: 


B = {'■ Nj :) — {Nj) 2 (j ~ s,i) 

where the symbol : : indicates normal ordering. In terms of the normalized g^ coefficient: 

(Nj ) 2 K 


(49) 


(50) 


In this way, as recognized in [221123] . the Schmidt number can be related to measurable statistical properties of light. 
In particular, formula (50) is well know to describe the statistics of multi-mode thermal light, with 1C playing the role 
of the ” degeneracy factor” characterizing the effective number of independent modes in a thermal beam. 

Figure [5] shows our results for the Schmidt number. The solid lines plot the ”exact” results, where 1C has been 
calculated by numerically performing the integrals involved in (46), (|47|), with the phase matching calculated via the 


complete Sellmeier relations. The red dashed lines in plot (a) and (b) are asymptotic behaviours, analytically derived 
by exploiting the linear approximation for phase matching. In particular, by using the approximated formula (35) for 
the coherence function, and performing the integrals involved in ( [46] ) and (47), one obtains the limit of the Schmidt 
number for a long pump pulse 


JC r '^T 


7T ^gvs 


2 V 2 AfL 


(51) 


For an ultrashort pump pulse, the asymptotic behaviour of 1C is calculated by using formula (39) or (37), for either 
the signal or the idler coherence function (identical results are indeed obtained). In this case 


K Tp ^T 


1 — T] V 7T Cl 


2 AIT 


fcgvm 


(52) 


The calculated asymptotes are well in accordance with our qualitative estimates of the number of modes in Sec |III[ 
based on the ratio between the spectral bandwidth and the coherence length. 

This shape of the curve, showing a minimum of 1C for a given value of the pump bandwidth and linear asymptotes 
at small and large values of the bandwidth, is commonplace, with a qualitatively similar curve characterizing also 
the co-propagating case in either temporal [5] or spatial [24] or even spatio-temporal [4 domains. The novelty here 
is that the minimum value of 1C is very close to unity, and remains close to unity for a rather large range of ACl p 
(see panel (c) in Fig(5]). This represents indeed a big difference compared to the copropagating case, where in order 
to generate separable biphotons very special matching conditions have to be chosen, corresponding to a zero group 
velocity mismatch between the pump and one of the twin photons, which can be realized only in type II interactions 


In the backward propagating case the conditions for separability are very easily approached, and rely entirely on 
the fact that rj = T gvm /Tg VS is naturally a very small quantity, because the temporal separationsr gvm , r gvs between the 
co-propagating and the counterpropagating waves are on well separated time scales. 

Indeed, a more refined calculation shows that the minimum value of /C, reached for a pump duration intermediate 
between r gvm and r gvs is K m i n = 1 + 0(g). This is also confirmed by analytical calculations of the Schmidt number, 
reported in detail elsewhere[27], performed by means of a Gaussian approximation of the sine function of phase 
matching, similarly to what done in [25] . These calculations (plotted as the blu dash-dot line in FigjSJ:) show that 
the minimum of 1C is 

)Cmin = ~ 1 + 2 V, (53) 

reached for Af 1 p = ^/3Cl' gvs Cl gvm This result suggest that a higher degree of purity of the reduced states can be 
achieved as the GVM between the two forward propagating is reduced. This is confirmed by the examples in Fig[6j 
which compares different crystals and phase matching conditions. Notice that a small GVM corresponds to a higher 
degree of purity, as in the copropagating case, but that in the present case the condition for separability is much less 
demanding, as it does not require a vanishing GVM, but just that that r gvm is small compared to sum of the inverse 
of group velocities Tg VS , which is always verified to some extent. 
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T p ( P S) 


Tp (PS) 


1.70 0.85 0.57 0.43 0.34 0.28 



T p (ps) 


FIG. 5: Schmidt number, as a function of the pump spectral bandwidth (lower axis) or duration (upper axis), (b) and (c) 
are insets of (a), showing the transition from high entanglement for a long pump r p Tg VS to an almost separable state for 
Tgvs r p Tgvm. The red dashed lines in (a) and (b) are the calculated asymptotic behaviours, the blu dash-dot line in 
(c) is the result of a Gaussian approximation. 4 mm PPKTP A in Fig[5J with Tg VS = 25.5 ps r gvm = 0.27ps, 77 = 0.01, other 
parameters as in Fig [3] 


V. INTERPRETATION: THE BIPHOTON CORRELATION IN THE TEMPORAL DOMAIN 


An alternative insight into the issue of separability vs entanglement is provided by the analysis of the biphoton 
correlation in the temporal domain. We consider 


= <ir(*»)4 ?Ut (ti)> = I^J 


d^i 

—f=e 

V2tt 


^i(Q s t s -\-Qiti) ik s (Q s )l c 


^(0,3 , r ^), 


(54) 


which is proportional to the probability amplitude of finding a signal and an idler photons at their crystal end 
faces at times t si ti. By using the linear approximation for phase matching (15) and performing the simple Fourier 
transformations involved in (54) we obtain: 


1 ti) — 


ge 


ik s l 


2 T, 


gvs 


sLc (- t s — (t s —ti 

&p [t s + r] — - ) Rect 


I-77 


2 TV 


gvs 


(55) 
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(ps) 

1.02 0.51 0.34 0.26 0.20 0.17 0.15 0.13 



FIG. 6: Role of GVM in determining the purity of the state: Schmidt number for different crystals and/or different phase 
matching conditions, corresponding to the points A, B, C in Fig[5] A) 4 mm KTP with r gvm =0.27ps, Tg VS =25.5ps, —>> rj = 0.01 
( same as in Figs§[5}.B ) 4 mm KTP, with r gvm =0, Tg VS =24.7ps, 77 = 0 . C) 4 mm LiNb 03 with T gV m=l-68ps, Tg VS =31.2ps, 
rj = 0.05 


where 


Rect(x) = 


1 


for 


xe 


0 elsewhere 


(-id) 


(56) 


is the box function of unitary width. The barred arguments t s ,ti denote time intervals measured starting from the 
arrival times of the centers of the signal/idler wavepackets. Precisely, t s ^ = t s ^ — tAs,i > where 


t As = (K + fcph at t Ap - (kp - k' s 
t A i = ( k- + k p )~ = t Ap - (k p - fc'F 


(57) 

(58) 


where tAp = k! p l c is the time when the center of the pump pulse exits the crystal slab. Figure [7] shows three examples 
of the temporal correlation function (55). 

The general formula ( [55] ) can be simp lified in the limit where the pump is long with respect to r gvm , i.e. in the 


quasi CW or intermediate limits ( 21 , 28k , where it takes the form 


,/T - x ' p ^ ' Gr V M -U 1 / — \ 

4>(ts,ti) — ge S c a p {t s ) 


Indeed, when the the pump pulse is long with respect to r g 


1 


2t £ 


-Rect 


gvs 


2t £ 


gvs 


(59) 


we have a p (t s + 77 t [_ t p ^ ~ &p(t s ) , because \t s — U 


is limited by the box function to values smaller than r gvs , so that 77 — t, 


— Tgvm pp 


Formula (60) shows that in the limit of a negligible GFM, the distribution of separations t s — ti between the 
arrival times of the twin photons is entirely described by the box function of width 2r gvs . This form of the temporal 
correlation clearly reflects the spontaneous character of the process, where photon pairs can be generated at any point 
of the crystal with uniform probability. Thus, assuming for simplicity that the the twins travel with the same group 
velocities v gs = v g i, the separation between their arrival times ranges with uniform probability from zero, when the 
two photons are generated at the center of the crystal up to ±r gvs = l c /v g , when they are generated at each of the 
end faces. 1 


1 Precisely, when the two photons are generated at the crystal center t s — ti = tAs ~ tAi = {k f s ~ ^' i )l c /2 « 0 , and the delay between 
their arrival times ranges uniformly between i) t s — ti = tAs ~ tAi ~ T gvs = — when they are generated at the right end face of the 
slab, and ii) t 3 — ti = tAs ~ tAi + T gvs = k' s l c when the photon pair is generated at the left end face. 
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(c) T p =0.06 ps 



FIG. 7: (Color online) Temporal correlation of twin photons \</>(ii,t a )\ , given by Eq.(|55|), plotted in the plane (a)High 

entanglement case, with JC ~ 26 , for a quasi CW pump r p = 14rg VS . (b)Almost separable case with JC ~ 1.06, for an 
intermediate pump r p = 0.04rg VS = 4r g vm. (c)Ultrashort pulse r p = 0.22r g vm, corresponding to an entangled state with JC ~ 4. 
Same KTP crystal as in Fig [5] 


The CW pump limit (Fig{7^t) corresponds to the situation where the pump pulse is much longer than the maximal 
temporal separation r gvs between the twins. In this case, the usual picture of the temporal entanglement of twin 
photons holds: the time when a signal or idler photon is individually detected has a large undeterminacy, because 
a photon pair can be generated at any time along the pump pulse. However, from the arrival time of one of the 
members of the pair one can infer the arrival time of the other with a much smaller uncertainty r gvs , which represents 
the mean uncertainty in the arrival time of one photon provided its twin have been detected, i.e. the correlation time. 
This kind of correlation is basically what predicted in Ref. [13] for a strictly monochromatic pump. 

However, when the pump pulse shorten below r gvs (Fig (7b) this description ceases to be valid, because the local¬ 
ization of the pump pulse provides a timing information on the arrival time of the signal that is more precise than 
the uncertainty in the temporal separation of the twins. Indeed when the pump pulse is much shorter than r gvs , but 
still long enough that GVM is negligible, the signal wavepacket overlaps almost exactly with the pump pulse during 
propagation, and the uncertainty in the arrival time of the signal is just the pulse duration. This is much smaller 
than the conditional uncentainty r gvs by which the arrival time of the idler can be inferred from that of the signal, so 
that the arrival times of the members of the pair appear completely uncorrelated. Indeed, the temporal correlation 
in Fig (TJ d is approximately: 


— 9^ S C( ^p (j's) 


2tv 


-Rect 


gvs 


2 Tc 


gvs 


(60) 


which is a factorable function of t s , U . 

Notice that when the pump pulse is so short that GVM starts to be important (Figj7j:), there is again a loss of 
absolute timing information, because in this case the arrival time of the signal cannot be inferred from that of the 
pump with a precision better than r gvm . In contrast, the arrival time of the signal conditioned to a photon count in 
the idler arm can be predicted within the short pump duration r p , and the state becomes again entangled. This can 
be better understood by looking the correlation function (55), which for r p T gvm can be rewritten as 


ti) 


ge ( ts - ijtj \ 

2 r gvs ap { I-77 J 
ge lkJc 

2r gvs V 1 - V ) 


Rect 


Rect 



(61) 

(62) 


where the last line has been obtained by substituting G — t s /v inside the argument of the box function (valid because 
the pump profile is much narrower than both r gvs and r gvm ). From formula (62) we see that, provided that an idler 
photon is detected, say at time G, the arrival time of the signal can be predicted as t s = rjti within the narrow 
uncertainty of the pump duration r v (see also Fig^). However when the idler is not detected, the overall uncertainty 
in the signal arrival time is the larger width r gvm of the box function. Clearly this argument predicts an entangled 
state, with the number of modes scaling as T gvm /r p , in agreement with formula (52). 
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VI. CONCLUSIONS 

In this work we provided a detailed theoretical analysis of the effect of the pump spectral properties on the quantum 
correlation of counterpropagating photons generated by SPDC in a periodically poled crystal. 

In particular, for increasing spectral bandwidths of the pump (descreasing pump durations), we demontrated a 
remarkable transition from a high-dimensional entangled state, to an almost separable state. The transition occurs 
when the pulse duration shorten below the characteristic transit time r gvs = ~ + 2 T - ^S^t along the 

crystal slab. This long temporal scale is a unique characteristic of the counterpropagating geometry, being associated 
to the delay between the times at which the counterpropagating photons, generated at some point along the slab, 
appear at their exit faces. The temporal correlation (temporal entanglement) is again restored for pump durations 
below the short temporal delay occurring between the co-propagating waves because of their different group velocities. 
The natural existence of such separated time scales ensures the possibility of generating high purity single photon (i.e. 
a separable two-photon state), under very general conditions, which differs drastically from the usual co-propagating 
geometry [26] . 

These conclusions have been supported throught the paper by the analysis of the Schmidt number in Sec |IV[ and 
by analytical and numerical evaluations of the spectral and temporal correlation function (SecJIT| SecJV|. 

The study of the maginal statistics of photons in Sec. m has revealed several non-trivial features: 

While for a long pump pulse twin photons have the same spectrum and the same coherence properties, in the regime 
of separability they exhibit very different features. In particular, the properties of the counterpropagating idler are 
entirely determined by the phase matching in the medium, so that we can say that they reflect the momentum 
conservation in the process. On the other side, the spectro-temporal properties of the signal are a replica of those of 
the co-propagating pump laser pulse, and rather reflect the energy conservation. 

For an ultrashort pump pulse, our quantum analysis has retrieved results analogue to what predicted in the classical 
description of the MOPO [HIE], but with some additional limitation. At difference with the MOPO prediction [H], 
our results impose a precise inferior limit to the observable bandwidth of the backward idler photon, which cannot 
be narrower than the phase matching bandwidth U gvs . Clearly, our analysis is limited to SPDC, but we notice that, 
to our knowledge, measurements of the spectrum of the backward wave in the MOPO have been limited by the 
spectrometer resolution m , so that our findings may open a question about the effective bandwidth of the backward 
wave. 


Appendix A: Approximations for the biphoton amplitude 


In this Appendix we derive the approximated forms of the biphoton amplitude used in the text, which holds in the 
various pump regimes. In all the cases we make use of the linear approximation for phase matching (15), based on the 


assumption that the bandwidths in play are narrow so that dispersion can be neglected. Under this approximation, 
the general espression (10) of the biphoton amplitude becomes 


= ~i=a p (U s + Lh)U 

v 27r 




n 


gvm 


Qi 

SI' 


gvs 


(Al) 


where for brevity of notation we introduced the phase matching function V(s) = sine ( s ) e ls . 

We consider first the limit of a CW pump ©. Since the pump bandwith is much narrower than the bandwidths 
Dg VS and U gvm of phase matching, the presence of the pump Fourier amplitude term forces = — Qi into the phase 
matching function. As a result 


lim 

T P / T gVS^°° 



a p (Q s + Sli) V 




OL p (Q s + Sli) V 



(A2) 

(A3) 


where we used the relation 1/U gvs = 1/Sl' — 1/U gvm , according to the definitions (16 -18). 


The limit (24) of an ultrashort pump is also straightforward. In this case the bandwidths of phase matching are 


assumed to be much narrower than the pump bandwidth Ug VS <C U gvm <C 


&p, so that the phase matching function 
has a narrow peak, which on the slow scale of variation of the pump forces = —rj£l s , or Q s = —Cti/rj inside the 
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pump argument. Therefore 


lim ip(Q s ,£li) = - 
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(A4) 

(A5) 


The intermediate pump limit (28) where to' gvs <C AQ P <C Q gV m is a bit more involved. We remind that the existence 
of this limit also requires r] = Hg^/f2 gvm <C 1, which is in practice always verified to some extent. By introducing the 
pump frequency to p = to s + , we recast the argument of the sine function 
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(A6) 


where the term to p /to gvm has been neglected because is on the order AQ p /Q gvm 
Concerning the pump amplitude we recast it as: 
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(A7) 


where the approximation in the second second line holds because ( —f- ) is the argument of the sine function 

V Si S vm ' ^gvs J 

(see Eq.(Al)), so that it is limited to values inside the bandwidth of the sine, say on the order 10 . Provided 
thar to' gvs /Ato p is small enough, this term becomes therefore negligible. With this in mind we can write the limiting 
behaviour of the biphoton amplitude: 


lim 'ipitog, toi) = _ 6l p (1 — 77 )] e ^s vm x sine 

T p/ r gvs^° V27T 

4gvm/ T~p ^0 
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(A8) 
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where the approximation in the last line is not mandatory, but could be useful in order to get consistent results, 
because clearly this limit can be realized only for 77 = r gVYn /r gvs —»• 0 . 


[1] C. K. Hong, Z. Y. Ou, and L. Mandel, Measurement of Subpicosecond Time Intervals between Two Photons by Interference, 
Phys. Rev. Lett. 59, 2044 (1987). 

[2] A. Gatti, E. Brambilla, L. Caspani, O. Jedrkiewicz and L. A. Lugiato, Phys.Rev. Lett. 102 , 223601 (2009). 

[3] O. Jedrkiewicz, J.-L. Blanchet, E. Brambilla, P. Di Trapani and A. Gatti, Phys. Rev. Lett. 108 ,253904 (2012) 

[4] A. Gatti, E. Brambilla, T. Corti and D. M. Horoshko, Dimensionality of the spatio-temporal entanglement of PDC photon 
pairs, Phys. Rev. A 86, 053803 (2012). 

[5] Y. M. Mikhailova,P. A. Volkov and M. V. Fedorov,Phys. Rev. A 78, 062327 (2008) . 

[6] Malte Avenhaus, Maria V. Chekhova, Leonid A. Krivitsky, Gerd Leuchs and Christine Silberhorn, Experimental verification 
of high spectral entanglement for pulsed waveguided spontaneous parametric down-conversion, Phys. Rev. A 79, 043836 
(2009) 

[7] S. E. Harris, Proposed backward wave oscillation in the infrared,Appl. Phys. Lett. 9, 114116 (1966). 

[8] Carlota Canalias and Valdas Pasiskevicius, Mirrorless optical parametric oscillators, Nature photonics 1, 459 (2007). 

[9] Gustav Strmqvist, Valdas Pasiskevicius, Carlota Canalias, Pierre Aschieri, Antonio Picozzi, and Carlos Montes, Temporal 
coherence in mirrorless optical parametric oscillators, J. Opt. Soc. Am. B 29, 1194 ( 2012). 

[10] Valdas Pasiskevicius , Gustav Strmqvist, Fredrik Laurell, Carlota Canalias, Quasi-phase matched nonlinear media: Progress 
towards nonlinear optical engineering, Optical Materials 34 513523 (2012). 

[11] Yujie J. Ding and Jacob B. Khurgin, Backward Optical Parametric Oscillators and Amplifiers, IEEE J. of Quantum 
Electronics 32 ,1574 (1996). 

[12] Tommaso Corti, Enrico Brambilla and Alessandra Gatti, Critical behaviour of the temporal correlation and coherence of 
counterpropagating twin beams, preprint. 

[13] Toshiaki Suhara and Makoto Ohno, Quantum Theory Analysis of Counterpropagating Twin Photon Generation by Para¬ 
metric Downconversion, IEEE J. OF QUANTUM ELECTRONICS 46, 1739 (2010) 





























17 


[14] A. Christ, A. Eckstein, P. J. Mosley and C. Silberhorn,Pure single photon generation by type-IPDC with backward-wave 
amplification, Opt. Expr. 17, 3441 (2009) 

[15] M. C. Booth, M. Atature, G. Di Giuseppe, A. V. Sergienko, B. E. A. Saleh, and M. C. Teich, Counter-propagating 
entangled photons from a waveguide with periodic nonlinearity, Phys. Rev. A 66, 02381502323 (2002). 

[16] E. Brambilla, O. Jedrkiewicz, P. Di Trapani, A. Gatti, Space-time coupling in upconversion of broadband down-converted 
light, J. Opt. Soc. Am. B 31, 1383 (2014). 

[17] Robert W. Boyd, Nonlinear Optics, 3rd Edition, Chapter 2, p.84 , Elsevier ( 2008) 

[18] David Nikogosyan, Nonlinear Optical Crystals: A Complete Survey, Springer New York (2005). 

[19] A. Ekert,P.L. Knight, American Journal of Physics, 63, 415-423 (1995) 

[20] S. Parker, S. Bose, and M. B. Plenio, Phys. Rev. A, 61, 032305 (2000). 

[21] M.P Van Exter, A. Aiello, S.S.R. Oemrawsingh, G. Nienhuis and J.P. Woerdman, Phys. Rev. A 74 012309 (2006). 

[22] K. Laiho, A. Christ, K. N. Cassemiro, and C. Silberhorn, Testing spectral filters as Gaussian quantum optical channels, 
Opt. Lett 36, 1476 (2011). 

[23] Andreas Christ, Kaisa Laiho, Andreas Eckstein, Katiscia N Cassemiro and Christine Silberhorn, Probing multimode 
squeezing with correlation functions, New J. of Phys. 13, 033027 (2011). 

[24] C. K. Law and J. H. Eberly, Analysis and Interpretation of High Transverse Entanglement in Optical Parametric Down 
Conversion, Phys. Rev. Lett. 92, 127903 (2004) . 

[25] Grice W P, U’Ren A B and Walmsley I A Phys. Rev. A 64 063815 (2001) 

[26] Mosley P J, Lundeen J S, Smith B J, Wasylczyk P, U’Ren A B, Silberhorn C and Walmsley I A, Heralded Generation 
of Ultrafast Single Photons in Pure Quantum States, Phys. Rev. Lett. 100 133601 (2008 ); P. J. Mosley, J. S. Lundeen, 
B. J. Smith, and I. A. Walmsley, Conditional preparation of single photons using parametric downconversion: a recipe for 
purity, New J. Phys. 10, 093011 (2008). 

[27] E. Brambilla, A.Gatti, The Schmidt number of counterpropagatin twin photons, in preparation 



